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We derive upper and lower bounds on the convergence behavior of certain classes of one-parameter quan- 
tum dynamical semigroups. The classes we consider consist of tensor product channels and of channels with 
commuting Liouvillians. We introduce the notion of Cutoff Phenomenon in the setting of quantum information 
theory, and show how it exemplifies the fact that the convergence of (quantum) stochastic processes is not solely 
governed by the spectral gap of the transition map. We apply the new methods to show that graph states can be 
prepared efficiently, albeit not in constant time, by dissipation, and give the exact scaling behavior of the time 
to stationarity. 



I. INTRODUCTION 

The Quantum Information and Computation community has thus far mostly focused on discrete time 
processes and unitary dynamics. Recently, however, the focus has been shifting towards more physically 
motivated processes, which inevitably involve continuous time and open systems dynamics. One of the 
major bottlenecks in understanding the information processing potential of open quantum systems has been 
the lack of quantitative tools for analyzing their convergence to stationarity. The interest in understanding 
the convergence behavior of quantum dissipative time-evolutions and quantum channels has cropped up in 
several distinct areas of quantum information theory, and seems to require a new set of tools for its analysis. 
The aim of this paper is to present one tool in this direction - the cutoff phenomenon - and illustrate it in 
explicit examples and applications. 

As a first step in developing tools for studying convergence rates of quantum processes, it is important to 
understand how much of the machinery from the classical Markov chain literature can be borrowed. Indeed, 
by noting that a finite-dimensional quantum channel is the non-commutative analogue of the probability 
transition matrix of a finite-state Markov chain, many results from the field of Markov chain mixing can 
be translated to the quantum setting, with appropriate modifications. A first effort in this direction was 
taken in iQl], where x^-convergence was considered, a quantum version of detailed balance was defined, 
and a restricted quantum Cheegers bound was proved. Furthermore, Hilbert's projective metric jld yields 
quantities which upper bound convergence rates. These studies and the few other attempts at quantifying 
the convergence behavior of quantum processes have focused on convergence for asymptotic times. 

More generally, it has become almost a folk theorem that "the gap governs the convergence" of (quantum) 
stochastic processes. This statement is true in the sense that the spectral gap represents the exponential 
convergence rate of a (quantum) stochastic process for sufficiently large times. However, there are situations 
where the asymptotic behavior either does not kick in before an amount of time exponential in the system 
size, or does not address the relevant physics. In particular, the important question is often "after how much 
time, as a function of some dimensional parameter (e.g. the system size), does the gap properly describe 
the convergence behavior?" In the case of state preparation |3], dissipative quantum computation [4], or 
thermalization |l5|], it is critical to guarantee that the asymptotic regime is reached in a time which is not 
exponential in the system size. In the case of quantum memories and error correction, on the other hand, 
one would like to guarantee the opposite direction, i.e. that the information initially encoded in the system 
is preserved for as long as possible, e.g. for an exponential amount of time. 

One of the few situations where one can make rigorous statements about the pre -asymptotic behavior of 
classical Markov chain convergence is known as the Cutoff Phenomenon Loosely speaking, the cutoff 
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phenomenon describes the situation where the (quantum) Markov chain, for some initial states, stays far 
away from its stationary distribution for a possibly long time (thus, e.g. preserving classical information), 
and then, at a specific time that may depend on the system size, suddenly approaches the fixed point (thereby 
suddenly losing all information from the initial state). Cutoff phenomena depend, apart from the type of 
Markov chains, on the chosen distance measure. The prototypical example of this behavior is in card 
shuffling 101, where a deck of cards is well shuffled only after a number of shuffles logarithmic in the deck 
size, whereas before that time it has large dependence on the initial ordering. In this article, we state a 
quantum version of this framework, and apply it to some situations of relevance in Quantum Information 
Theory. 

The outline of the paper is as follows. In Section|IIl we set the notation and introduce basic mathematical 
tools, in particular spectral properties of quantum dissipative time evolutions and distance measures. In 
Section|IIIl we give a definition of the Cutoff Phenomenon in the quantum setting. In Section|IV] we state 
and prove the main results, namely cutoff-type bounds for time-evolutions due to commuting Liouvillians 
or tensor product channels. Section[V]illustrates these main results by various examples. As an application, 
we show in particular that the dissipative preparation of graph states takes a time logarithmic in the system 
size, and is thus not solely governed by the spectral gap of the Liouvillian. Finally, the conclusion and a 
brief outlook are given in SectionfVll 

II. BASIC NOTIONS 
A. Mathematical Setting 

Throughout the paper, we restrict ourselves to finite-dimensional quantum systems. The set of 
quantum states (density matrices) acting on a d-dimensional Hilbert space is denoted by Sd = 
{p e -Mdlp ^ p\p> 0,tr[p] = l}, its restriction to full-rank density matrices by . Quantum chan- 
nels are completely positive trace-preserving linear maps T : Md ^ where here the input and output 
spaces are taken equal to allow for repeated application of the channel. The dual map, i.e. the channel in the 
Heisenberg picture, is denoted T* and by definition satisfies tr[^r* {B)] = tv[T{A)B] for all A, B e Md- 
We consider almost exclusively quantum channels Tt :~ e*^ (for t > 0) that are elements of the quantum 
Markov semigroup generated by a Liouvillian C, which due to the trace-preserving property of Tt satisfies 
C*{1) = 0. Tt describes the evolution after time t coming from the master equation p = C{p). When fixing 
some complete orthonormal basis {-F'j}i=i,...,ci2 of Md, any linear map S : Md Md can be represented 
as a X d^-matrix S by defining Sij := tv[F} S{Fj)]. We dress the operator S with a hat when working 
in the matrix representation. 

B. Spectral Properties of Quantum Channels 

The spectrum S(r) of a quantum channel T : Md A^d is contained in the closed unit disk of the 
complex plane, and every eigenvalue has a complex conjugate partner One eigenvalue is 1, the correspond- 
ing eigenspace is spanned by quantum states, and the density matrices in this eigenspace are exactly the 
stationary states of the channel. Like any linear map, a quantum channel can be decomposed uniquely into 
Jordan normal form: 

f = ^ PkPk + Nk , 

k 

where each fik G S(T) is an eigenvalue of T, Pk is a projection onto the corresponding Jordan eigenspace, 
and Nk is a nilpotent matrix on this eigenspace. It can be shown that the eigenvalues pk with modulus 
\pk\ = 1 (the so-called peripheral spectrum) have trivial Jordan blocks {Nk = 0), and thus the asymptotic 
space of the quantum channel T is Ht := span{X G A^^lEI^ £ IR : T{X) = e'^X}. Its asymptotic 
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evolution is the phase-preserving projection onto H^, 

k:\fik\ = l 

and the corresponding map T^p : A4d — > is a quantum channel. If Tt = e*^ comes from a one- 
parameter semigroup of quantum channels, the spectra and (generalized) eigendecompositions of Tt and of 
C are related by exponentiation. In this case, we suppress the time-dependence when writing := {Tt)ip. 

If a channel T has only one eigenvalue of modulus 1 and if its unique stationary state has full rank, then 
we call T primitive We call a Liouvillian generator C primitive iff e*^ is a primitive channel for all 
t > or, equivalently, iff e*^ is primitive for one t > or, equivalently, iff C has exactly one eigenvalue 
and the corresponding stationary state p (i.e. C{p) = 0) has full rank. 

Note that [T — T^) has spectral radius equal to jl where jl sup{|A| | A € S(r), |A| < 1} is the largest 
modulus of the eigenvalues of T in the interior of the unit disk. If Tt — e*^ comes from a one-parameter 
semigroup, then /l(t) = e~*^ where A := inf{|Rc(A)| | Re(A) < 0, A e !](£)}, and A is referred to as 
the sap of the Liouvillian. For a more detailed discussion of spectral properties of quantum channels, see 

C. Distance Measures and Contraction Coefficients 

There are many ways to quantify distance between two quantum states. We recall basic relationships 
between the distance measures used in this paper. For p,<T G Sd define the distance measures: 

1. The trace distance: (itr(p,o') = |||/9 — crU i, where ||X||i = tr[vXTX] is the trace-norm. 

2. The Bures distance: dsip, cr) = ■\/l — F{p, a), where F{p, a) = ^'^'[\/^/p^\/p] is the fidelity. 

The trace and Bures distances have the important property that they are monotone under the application 
of quantum channels, i.e. d{T{p), T{a)) < d{p, a) for any channel T and any states p, a. It is well known 
that they satisfy ifTTIl : 

dl{p,(T) < dtr{p,a) < dB{p,a)^2-dl{p,a) = ^1-F^p,a) . (2) 

Moreover, for states p close to a full-rank state a, the Bures and trace distances can be bounded even linearly 
in terms of each other (see Appendix lAlfor a proof): 

Proposition 1 (Trace vs. Bures distance infinitesimally) Let a € be a strictly positive density opera- 
tor with smallest eigenvalue Xminicr) > 0. Then there exists e = e(a) > such that 

dB{p,cr) < -j=l===dtr{p,a) (3) 

for all density matrices p G Sd with dti{p, a) < e. 

Having introduced appropriate distance measures between quantum states, we now want to characterize 
the distance to stationarity of a quantum channel. Such a convergence measure should somehow quantify 
one or several of the three closely related properties: (i) how far an output state of the channel T may be 
from its limiting evolution T^p, (ii) how reversible the action of the channel is, (Hi) how much information is 
lost by the application of the channel. We choose the following definition, which addresses point (i) above. 

Definition 2 (Trace-norm contraction) Let T : Aid ^ -Md be a quantum channel. Then its (trace-norm) 
contraction is defined as 

r,tr[T] := sup dtr (r(p), T^(p)) . (4) 
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The Bures contraction 775 [T] is defined analogously by using the Bures distance ds- Eqn. (01) is not the only 
possible choice of a contraction measure, and like the other possibilities it has strengths and weaknesses. 
We point out in particular that Eqn. ^ can be discontinuous in T, whenever the size of the peripheral 
spectrum is discontinuous. 

In the setting of dissipative state preparation (where the peripheral spectrum is usually engineered to 
be trivial), 77tr[r] equals the maximum trace distance between the desired steady state and the channel 
output, i.e. for the most disadvantageous^ chosen initial state p. Having this in mind, we now state an 
asymptotic convergence theorem for one-parameter semigroups of quantum channels, which is also a crucial 
fact guiding our further investigations. 

Theorems (Contraction theorem) Let C : Aid ^ M.d be a LiouvilUan, i.e. the generator of a one- 
parameter semigroup of quantum channels Tt = e^'~ [t > 0), and let A be the gap of C. Then, there exists 
L > and for any < X there exists R > such that 

Le-^^ < T]t,[Tt] < Re-^" Vt > . (5) 

A proof is supplied in AppendixiBl showing in particular that one may not choose ly = X when an eigenvalue 
Afe of £ with modulus |Afe| = A has a non-trivial Jordan block. Theorem[3] makes a statement only about 
the asymptotic convergence behavior, i.e. the exponential rate as t — s- 00. If is taken close to A, then 

R can become very large. But for fixed i^, a universal dimension-dependent upper bound on R of order 

,2 

d can be obtained. If the system at hand is composed of many particles, say n qubits, then this time- 
independent prefactor can in principle become doubly exponentially large in the number of particles; this 
would correspond to an exponentially long time to convergence, even for a gap constant in the system size. 

Finally, we note that an analogue of Theorem[3]holds also for the repeated application of a discrete-time 
channel T, and results analogous to the ones following can be formulated in the discrete-time setting as 
well. 



III. THE CUTOFF PHENOMENON 

Often times, the relevant question when considering the convergence behavior of open systems is "how 
long does the process have to run before it reaches equilibrium?" To make precise statements about con- 
vergence, it is usually necessary to consider how the time to convergence scales with the system size. For 
instance, one would like to know how fast, as a function of the lattice size, a given dynamical process on a 
lattice converges to its steady state. 

It was observed a while ago in the setting of classical Markov chains, that for a special set of chains this 
question can be answered exactly when the size of the system becomes large. This behavior, which has 
been coined the Cutoff Phenomenon [6-8], characterizes the situation when for some (possibly long) period 
of time some information from the initial state is perfectly preserved until a critical time. Shortly after this 
critical time, however, essentially no information of the initial state can be recovered from the time-evolved 
state anymore. For large system sizes n, the contraction of the channel as a function of time t will look like 
a step function at the cutoff time tn (see Fig.[T^). 

This behavior has been observed and proved to occur in a number of interesting examples of classi- 
cal Markov chains. One case where this phenomenon is particularly pronounced, and which triggered 
widespread popular interest, is in card shuffling, where it was shown that a deck of n cards is well mixed 
after exactly 3/2 log 71 riffle shuffles, and poorly mixed under 3/2 log n riffle shuffles (when n becomes 
large) In particular, this guarantees casino owners that if their dealers riffle shuffle their 52 poker cards 
seven or more times before each draw, then they do not need to worry about about players trying to improve 
their odds by counting cards'. Several other processes have also been shown to exhibit cutoffs, including 



The punchline in this example is that the cutoff behavior can be exploited to perform a spectacular magic trick called "Premo" fl^ . 
where the magician can guess a card chosen randomly by a person from the audience who afterwards shuffles the deck a few times. 
The infomiation in the initial configuration (i.e., the identity of the card that the person put on top) is preserved and can be identified 
before a "time" of exactly ^ log n riffle shuffles (for large n). 
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FIG. 1: |(a)| Behavior of the contraction r;tr[r/"'] as a function of evolution time t for a one-parameter semigroup of 
channels that exhibits a cutoff (for large system size n). At time t — t^, essentially all dependence on the initial state 
vanishes in a small window of time o(t„ ). |(b)| In the case of pre-cutoff, the contraction will be close to 1 for times 
i < = 6(fcn) and close to for t > t2,n = 0(fc„), but now the gap between ti^n and t2,n may be of order 0{kn). 



random walks on graphs with a discrete group structure, birth and death type chains, and some Monte Carlo 
sampling methods Isl [Tsi [l4l l . We now give a formal definition in the quantum setting: 

(n) 

Definition 4 (Cutoff) Let ' be a sequence, indexed by the "system size " n, of one-parameter semi- 

(n) 

groups of quantum channels. We say that exhibits a cutoff (in trace-norm) at times tn, if for any real 
c > 0: 

c < 1 => lim ?7tr[ri"j] = 1 , 

n— >cx3 

Ol ^ lim r;tr[rir'] = . 

We point out that this is not the only definition of a cutoff, and in a sense it is an incomplete one, as it does 
not provide detailed information about the cutoff window, i.e. the width of the drop-off at t„ (see Fig.[T^). 

In many situations, it is difficult to prove an actual cutoff, whereas it might be easier to show a weaker 
statement which gives only the precise order of magnitude of the time to convergence: 

Definition 5 (Pre-cutoff) With the same assumptions as in Definition we say that T^^^ exhibits a pre- 
cutoff (in trace-norm) of order Q(kn) for some sequence kn, if there exist times ti^„ and t2,n that are both 
of order 0(fc„), such that 

c<i ^ lim ,7t,[Ti;;>j 1 , 

Ol => lim Tjt,[Tj:^^ ] = 0. 

We again emphasize that the cutoff phenomenon does not cover all types of pre-asymptotic behavior 
For example, a process can follow a polynomial decay for a certain amount of time before it settles into the 
asymptotic (exponential) regime. We also point out that the Cutoff Phenomenon can be defined with respect 
to any monotone distance measure (even unbounded ones like the x^-divergence). A cutoff in one distance 
measure does not imply a cutoff in another measure. This suggests that obtaining such a tight contraction 
estimate might in some cases reveal information about only one specific facet of the convergence behavior, 
associated with the given distance measure. For instance, the trace-norm contraction measure of a channel 
tells us whether a single bit of classical information is preserved after a certain time, but makes no statement 
about the entire amount of information or the preservation of quantum information. 
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rV. MAIN RESULTS 

In this section, we present two situations which exhibit behavior related to the Cutoff Phenomenon intro- 
duced above: commuting Liouvillians and, as a more restrictive but still relevant class, Liouvillians acting 
independently on subsystems. In the latter situation, the contraction behavior of the constituent channels 
is known, and we ask how the contraction of the tensor product of channels behaves. We show in general 
terms that a sequence of tensor product channels exhibits a pre-cutoff of order 8(logn), where n is the 
number of tensor factors. In the next section, we discuss some specific situations of this kind where an 
actual cutoff occurs, which includes the dissipative preparation of stabilizer states. 

Our first theorem provides a general upper bound on the contraction of a channel from a one-parameter 
semigroup whose Liouvillian is composed of commuting parts, i.e. C = J^j where [Cj , Ck] = 0. In this 
case, the gap of the full Liouvillian is at least the minimum of the gaps of its constituent parts. We show 
that in this context the time to convergence is upper bounded by O(logn) times the convergence time of 
the "slowest" constituent channel, where n is the number of commuting terms in the Liouvillian. 

Theorem 6 (Contraction for commuting Liouvillians) Let Cj : Md ~^ Mdbe Liouvillians which com- 
mute, i.e. [£j, Ck] = 0/or all j,k = 1, . . . , n. Define C = Cj, and the corresponding semigroups of 

channels Ttj = e^^' andTt = e*^ [t > 0). Then: 

ritATt] < E'?*' !^*.^] • (6) 



Proof The theorem is proved by induction. Let T^pj be the projector onto the asymptotic space of Tt i, 
and let T^j^i be the projector onto the asymptotic space of Ttj^i = e*^3^i , see Eqn. ([T]i. Note that 
Tt = Tt^iTtj^i and = T^^iT^.j^i, by commutativity of the Liouvilhans. Then, 

VtATt] = ^ sup IKTt^iTt.j^i -Ty,iT^j^i)(p)||i 

= t; sup ||Tt_i(rtj^i - T^,j^i){p) + {Tt^i - T^s){T^,j^i){p)\\i 

< Vti[TtA] + Vtr[Tt,j^i] , (7) 

where the last inequality follows from the triangle inequality, from monotonicity of the trace-norm under 
quantum channels, and by definition of ?7tr [Tt] . By induction, we get r/tr [Tt] < J2j Vti [Tt,j]- □ 
An immediate consequence of Theorem |6] also using Theorem [3] is that for a system described by n 
commuting Liouvillians with bounded gaps, the convergence time will be upper bounded by O(logn). 
Note however, that commuting Liouvillians should not be confused with classical processes. Indeed, as 
described in Section IV Dl graph states can be prepared dissipatively as stationary states of commuting 
Liouvillians, whereas these states can be highly entangled (e.g. cluster state). 

Our second result gives the general contraction behavior of a tensor power of a one-parameter semigroup 
of quantum channels. This is a special case of commuting Liouvillians, but where strict upper and lower 
bounds can be derived as the number of tensor factors becomes large, thereby estabhshing a pre-cutoff as 
defined in Section HiH 

Theorem 7 (Pre-cutoff for tensor powers) Let C : Md Md be a Liouvillian with gap A, and let 

Tt = e*'^ [t > 0). The sequence of one-parameter semigroups T^'"'' = T®" exhibits a pre-cutoff in 
trace-norm at times ii,„ = log (n) /2A and t2.n = log (n) / A. 

Proof Here and below we use the fact that [T ® S)^ ^ T^p ® S^p for any pair of quantum channels T and 
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S. To prove the lower bound, let c € (0,1): 



sup dtr(riri(p),r("Hp)) > supdtr((T,,,,„(a))«",(T^(a))«") 



> 1 — cxp 

> 1 — cxp 



l-n sup dl (Tcti „(cr),T^(cr)) 



ne 

2 



-2cti „ A 



1 — cxp 



2 



1 (n — > oo) 



(8) 



The first inequality is obtained by restricting the supremum to product states p ~ a®^, the next from 
Lemma [8] (see below), and the last follows from Theorem |3] (with some constant L > 0). Hence, 



lim„^oo VtrlT^tL] = 1^ for c e (0, 1). 



For the upper bound, we apply Theorem |6] to get ryti [Tj ] < nrytr[7t ® idn-i], where id„_i is the 
identity channel on n — 1 sites. In the following paragraph we show jytriT* ^ idn-i] < 4:dr]tr[Tt]. Now, for 
any given c > 1 one can choose < A such that ciy/X > 1, and by Theorem[3]one can find R such that 
Vtr[Tt] < Re^"*^ for all t>Q. Combining all this, we finally get that for any c > 1, 



(n) 

Ct2,„ 



< AdnRe 



ARdn^-""/'^ ^0 (n ^ oo) 



(9) 



It remains to show that yytr [T id] < Ad'qt^ \T] for any channel T : ^ Md and any identity channel 
id : Md' — > M.d'- The inequality (used below) between the norm || • ||i_i on superoperators induced 
by the trace-norm and its stabilized version || • ||c;, is proven in jTsil (exercise 3.11). In the following, 
X g M.d ® -Md' and A + iB £ Aid denote arbitrary matrices, A and B are Hermitian, and P,Q G Aid 
are positive semidefinite with PQ = 0. Further note that ||A||i = \\{A + iB) + {A - iB)\\i/2 < 
i\\A + iB\\i + \\A-iB\\i)/2 = \\A + iB\\i, and similarly ||B||i < \\A + iB\\i. Thus: 

77tr[r®id] = sup dtAT<E)id{p),T^(E)id{p)) < sup i||((r - T^) ® id)(X)||i 



< ^\\T-T^\U < ^||T-r^||i_i = ^ sup \\{T-T^){A + tB)\U 



ll-Y||i<i 
d 
2 



||A+iB||i<l 

<d sup \\{T^T^){A)\\i = d sup ||(r-r^)(p-g)||i 

||A||i<l l|P-Qlli<i 

< 2d sup ||(r-r^)(p)||i =Mtju-[t]. 
\\p\\l<i 

□ 

Theorem|2]can be generalized to certain cases where Tj*^"^ is the tensor product of a set of one-parameter 



semigroups Tj that are not all identical. See 111611 for the analogous classical result. 

Theorem|2]establishes pre-cutoff rather than actual cutoff behavior But at the end of subsection lV Cl we 
show examples where, for any chosen ?- e [1, 2], a cutoff occurs at times t„ = log(n)/rA. This means that 
ti „ and t2,n in Theorem|7] viewed as upper and lower bounds on the contraction, are tight when expressed 
in terms of the gap A. 

The following Lemma completes the proof of Theorem |7] (cf. Eqn. (O); the inequalities ( fTOl i for the 
Bures distance will be used to show cutoff in Proposition|9] For collections pi , ai E Sd of density matrices 
(i = 1, n), define p^"' = (^"^^ Pi and similarly (t'-"\ 

Lemma 8 (Distances between tensor product states) Let Pi,ai € Sd, i = 1, ...,n, and denote by dtr, ds 
the trace and Bures distances, respectively. Then the following inequalities hold: 



1 — exp 



^(i|(pi,cri) 



i=l 



< dUp^^\a(-^) < 



^d|(pi,cr,) 



(10) 
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< dtr(p("\a(")) < ^dtr(p.,aO . (11) 

i=l 

Proof The fidelity is multiplicative under tensor products, cr'"') — n"=i ^iPi^ ^i)- Also, by 

induction it is easily seen that (1 — Xi) < ~ ^0 ™y collection of reals Xi G [0, 1]. Thus: 

n n n 

z— 1 i—1 i—1 

Since Y\ - e^^~^ > Y\ - xi whenever xi > 0, we get the lower bound in < fTOb : 

n n 
1=1 i=l 

The upper bound in ( fTTT i follows by a calculation similar to the one yielding Eqn. while the lower 
bound follows from ( fTOl i and two of the inequalities in (|2]l. □ 

V. EXAMPLES OF THE CUTOFF PHENOMENON AND APPLICATIONS 

Theorem |7] establishes a pre-cutoff and thereby estimates, up to a factor of 2, the time to convergence. 
The next natural question is: when does an actual cutoff occur? We discuss two such situations. The first 
concerns tensor powers of primitive channels where the input states are restricted to be separable, the second 
concerns tensor powers of channels whose unique fixed point is a pure state. In subsection C we give an 
explicit example of this situation (the qubit amplitude damping channel), and in subsection D we apply this 
to the dissipative preparation of graph states and discuss how the cutoff phenomenon determines the exact 
convergence time of this process. 

A. Primitive Channels with Separable Initial States 

Beyond Theorem |7] we can establish a sharp cutoff for primitive Liouvillians when the inputs are re- 
stricted to be fully separable quantum states between the n channels: 

■= { T.P^P' ^ ■ • ■ ^3 Pn I Pfc > 0, ^Pfe - 1, e } . (12) 

k k 

Proposition 9 (Primitive Liouvillians with separable inputs) Let C : Md ^ Mdbe the generator, with 
gap A, of a one-parameter semigroup of primitive channels Tt = e*^ {t > 0), and define the trace-norm 

contraction ofT^"''^ = T®''^ restricted to separable input states: 

CP[ri")] := sup dtr(Tr(p),T«"(p)) . (13) 

Then, the sequence of one-parameter semigroups T^""^ = T®^ exhibits a cutoff ( with respect to the con- 
traction measure rf^"^) at times t„ = log (n) /2A. 

Proof From the primitivity of the channel we get that Tip{p) = cr for any input state p G Sd, where a is 
the unique stationary state of C. Further, as a is of full rank, the new bound from Proposition[T]and bounds 
from Eqn. (|2j together with Theorem [3] show that, for any i/ < X, there exist constants R > L > such 
that Le"*^ < 775 [T*] < Re-*"" for all t > 0. 

The rest of the proof follows the same lines as the proof of Theorem |7] However, we first show the 
theorem here for the Bures metric, i.e. replacing dti in by ds. For proving the lower bound, the same 



1 — exp 



1 " 

-n'^dl{p^,a^) 



n 

1=1 
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arguments as the ones leading to Eqn. ^ show that lim„^.oo iJ^B^iTct^] = 1 for c G (0, 1). For the upper 
bound, note that due to convexity of the Bures distance (derived from concavity of the fidehty ifTHl ) the 
supremum in (T3[ is reached for a product state p = <S)^i=i Pi^ so that Lemma|9]can be appUed: 

(vTiT^:!]) = sup 4((g)ret„(pO,aM < 5] sup 4(7^ct„(p,),^) 

As in the proof of Theorem |7] for each given c > 1 one can choose ly and R accordingly to show 
lim„_j.oo VB^[Tct„] = 0, which proves a cutoff at times t„. Finally, by Eqn. Q, a cutoff in the Bures 
contraction rjg^ is equivalent to a cutoff in the trace-norm contraction ril°^, at the same times i„. □ 

We do not know whether the separable input assumption is actually necessary for Proposition |9] If the 
assumption were indeed necessary, then the statement would imply the possibility of increased storage time 
of classical information due to an entangled encoding. 



B. Channels with Unique Pure State Fixed Point 

Proposition 10 (Unique pure state fixed point) Suppose that the pure state tjj = |^')(V'I S is the 

unique stationary state of the Liouvillian C, which has gap A and generates the channels Tt := e^^ {t > 0). 
Then, Tj'"-* = T^**" exhibits a trace-norm cutoff at times in ~ log (n) /V, for some \ < D < 2A. 

Proof Since C has only one stationary state, the peripheral spectrum of e*^ is trivial for all i > 0, so 
that T^{p) = ijj for all p £ Sd- This, together with well-known inequalities relating the fidelity and the 
trace-norm between a pure and a mixed state d, yields: 



1- inf F\Tf"'{p),ip'^'') < ritr[Tf"'] < Jl- inf F2(rf"(p), V®") ■ (14) 
The last infimum can be evaluated explicitly in the case at hand: 



inf F2(rr(p),V^") = inf tr[Tf"(p)V'^"] = inf tr[p (T^V-)) . 
peSd" peSart peSdn 

= A„,„, ((r;(^))^") = [A„„;„(r;(v.))]" 

= [i-A_,(i-r;(v))r = (i-||r;(i-0)iu)", 

where in the last step we used that Tt is trace-preserving (rj*(l) = 1) and that Tj* is a positive map. Thus, 
from equation (fT4l) above: 



i-(i-||t;(i-v)||oo)" < vtATfl < v/i-(i-im*(i-V')IU)". (i5) 

Since T^p{p) ~ ijj for all p G Sd, we have T^{A) = ijjir[A] for all A G Aid, &o that its dual is given by 
T*{B) = ltr[Bi/)] for B e Md- Thus: 

||t;(i-v^)IU = ||(r;-T;)(v;)||oo. 

Now we consider this last equation in a similar way as in the proof of Theorem [3] in Appendix IB] When 
one writes ip as a linear combination of the generalized eigenvectors of £*, then considering large times t 
will essentially pick out the Jordan-eigenvalue(s) occurring in ip which has largest real part — P < (i.e., 
not the eigenvalue 0), and among these it will pick out the polynomial(s) of highest degree J > that are 
occupied (i.e., occur with non-zero coefficient in the linear decomposition of -ip)- This means that, for any 
arbitrarily chosen to > 0, there exist constants < Ci < C2 such that 

C.iPtye-''' < ||(t; - T*m\\^ < C2{vtYe-'' Wt > to . 
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Using this in Eqn. (flSl l gives 

This proves a cutoff for Tf'" at times t„ = (log n) / D, since for any constants c,K>0 and J > 0: 

lim l-(l-/^(:>ci.)^e— = 1- lim fl + '"^'^^^"g"^""'"^ 

= 1— lim exp (— _K'(clogn)'^n"'^~'^) = 

■n— >-oo 

As (— i^) is the real part of an eigenvalue of £, it is evident that P > A, and Theorem|7]shows v < 2A. □ 

Both infima in the upper bound and in the lower bound in Eqn. ( fT4b are attained for some pure product 
state p ~ (p®"^, even though this is not clear for the supremum that achieves rjtr^f"']- The paradigmatic 
example in the following subsection saturates the upper bound v = 2A, but we also provide modifications 
of this example where v takes on any values between A and 2A. 



r 1, c< 1 

\ 0, c> 1. 



C. Qubit Amplitude Damping 

The amplitude damping process (on qubits) describes the situation where the excited state |1) decays 
into the ground state |?/') := |0) at a constant rate 7. This corresponds to a Master equation with a single 
Lindblad operator L := .y/7|0) (1| and no coherent contribution: 

C{p) := LpLt - Il^Lp-\pL^L = 7 (^|0)(0| • (l|p|l) - - . (16) 

A straightforward calculation shows that: 

, tc^_ / e"''* , < t < (log2)/7 (attained for p = |1)) 

J~ \ e-''*/V^4(l - e--^*) ,i>(log2)/7 (attained for p oc |1) + Vl - 2e-r*|0)) . 

Thus, ritr[e*''''] decays asymptotically in time as e and not as e"''*, which one would expect from 

the analogous classical noise process (the classical Markov map has eigenvalues and —7, whereas the 
Liouvillian (fT&t has two additional eigenvalues —7/2). 

By Proposition[TO] the semigroups (e*-'-)®" exhibit a trace-norm cutoff at times t„ = (log n)/v for some 
V with A = 7/2 < i> < 2A = 7. V can be computed explicitly by using ||T(*(1 — "0)1100 = e in 
Eqn. ( fTSb . or from the proof of PropositionfTOlbv writing i/j as a linear combination of eigenvectors of £* : 

^ = |0)(0| = (17) 

where 1 and |1) (1| are eigenvectors of C* with eigenvalues and —7, respectively. Thus, = 7 = 2A, and 
the cutoff occurs at times t„ = (log n) 

For system size n it thus takes time 0(log n) before convergence happens, even though the Liouvillians 

£(") = /:(K)id«)...®id + id(g)£(8)...(8)id+ ... + id® id®. ..«)£, (18) 

which generate the semigroups T^^^ = e*''''"', have a gap A^"^ = A = 7/2 which is independent of n. 
Therefore, this example refutes the conventional wisdom whereby "the gap governs the convergence time". 

If, in addition to the Liouvillian ( fT6] l, there are also processes with Lindblad operators -y/a|0)(0| and 
V^|l)(l| acting on each qubit, then the steady state %}} and its decomposition (fTTI i into eigenvectors of the 
dual evolution operator are as above (in particular, v = 7), and a cutoff occurs at times i„ = (log n) j^. In 
this new situation, however, the gap is given by A = min{7, (7 + a + /3)/2}, which shows that the bounds 
on the cutoff time given by PropositionfTOl and implied by Theorem|2]are tight. 
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D. Dissipative Preparation of Graph States 

As an application of the results from subsections B and C, we consider the dissipative preparation of 
graph states. This task was considered in IHI^l, where it was shown that a set of local Lindblad operators 
can be constructed in a way analogous to Eqn. (fTST i, such that the unique stationary state of the process is 
the desired graph state and that the spectral gap of the process is independent of the number n of particles 
or stabilizer operators. 

We complete this analysis by showing that the convergence time, as measured by the trace-norm contrac- 
tion, scales as log n with the system size. The trace-norm contraction is the relevant quantity to consider in 
this case, as it quantifies the maximal failure probability when the graph state is used for further quantum 
information processing, like the cluster state for measurement-based quantum computing. We actually show 
that the dissipative preparation of graph states exhibits a cutoff (in trace-norm) at times of order 0(log n) in 
the sense of Definition |4] Although still efficient, the logn scaling of the preparation time again indicates 
that the gap does not fully determine the convergence behavior. 

Proposition 11 (Dissipative preparation of graph states) Any graph state on n sites, associated to a 
graph of maximal degree k, can be prepared dissipatively in a time of order log n by using n Lindblad 
operators that are at most k-local. 

In fact, for large n, the preparation procedure described in fi^ 0/ takes exactly time t„ ~ (log n) /j to 
converge to the desired graph state, where 7 is the decay rate ("f/2 = spectral gap) of each local Lindblad 
operator and when starting from the most disadvantageously chosen initial state. 

Proof Given a set {Sk}k=i of stabilizer operators, the unique state which is an eigenstate of Sk with 
eigenvalue +1 for every k is called a stabilizer state. Graph states UtIi are a special case of these and can 
be described by an undirected graph with n vertices. The stabilizer operators of the graph state are then 
Sk = njenbhd(fc) '^j' where nbhd(/c) denotes the set of all vertices connected to vertex k by an edge. 

The stabilizer operators of a graph state uniquely define a "graph basis", written as {|^ii,...,i„)}i,e{o,i}' 

byS'fc|$ji t„) = These basis vectors satisfy (Tf,|$jj^...^j^=i,...,i,J = |$ii,...,ifc=o,...,j„>, 

and the "graph state" is |$o o)- 

Define the n Lindblad operators IH (/c = 1 , . . . , n) 

Lk = Vl<^k^, (19) 

and observe that Lfc|$ji,...,j^=i ij = ^/J\^^l,...,^^,=o,...,^J and I/fc|^ii,... ,4^=0,... = 0. Thus, in the 

graph basis, each of these Lindblad operators acts as one term of the sum ( fTSl l acts in the computational 
basis. Therefore, together they act like the tensor product of amplitude damping channels in subsection C, 
now with the graph state as the stationary state. Proposition[TO]or, more explicitly, subsection C thus prove 
a cutoff at times (log n)/^ for the preparation of graph states. □ 

Note in particular, PropositionfTTIshows that, for the procedure described by Eqn. (T9[ . there exist some 
initial states for which one can guarantee convergence not to occur before time (log n)/'-f. 

VI. CONCLUSION 

In this paper, we have introduced the notion of the Cutoff Phenomenon in the context of quantum infor- 
mation theory and applied it to analyze the convergence behavior of some composite quantum processes 
in continuous time. In particular, we show that the convergence, measured in the trace-norm, of a tensor 
product of one-parameter semigroups of time evolutions always exhibits cutoff-type behavior We identify 
two specific cases (primitive channels with separable initial states, and channels with a unique pure fixed 
point), which exhibit a true cutoff. We provided a full analysis of the problem of dissipatively preparing 
graph states, and show that the convergence time scales as log n in the number of stabilizer elements ??. 

Finally, we conclude by noting two directions where the methods introduced in this paper could be of 
use. The first is the task of passive error protection in the presence of local noise. It was shown recently 
ifisll that, if the noise is locally depolarizing and allowing for arbitrary Hamiltonian control, an optimal 
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protection time of order O(logn) can be achieved. Theorem |7] gives a strict upper bound on the amount 
of time that one bit of classical (and hence also quantum) information can be encoded into n qubits, when 
every qubit is subjected to local noise, and no Hamiltonian control is allowed for. The upper bound happens 
to coincide with the one in |18!], indicating that their result might not be restricted to depolarizing channels, 
but could be a general feature of tensor product channels. Along similar lines, a second extension of the 
above results is in the study of continuous time quantum information theory, where channels are replaced by 
one-parameter semigroups, and standard objects, such as channel capacities and compression rates, become 
functions of time. 
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Appendix A: Proof of Proposition [T] 

We now give a proof of Proposition[T]from Section HTCl 
Proof Denoting by and \i) the eigenvalues and -vectors of cr (i = 1, d), it is shown in ifTgl] that 
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using that all norms on a finite-dimensional vector space are equivalent. The last term is upper bounded by 
Kdlj.{p, a) with some constant K = K{(t) > 0, and we bound the denominator by 2A,„i„(cr): 

rfB(p,cr) < T 2^ 77\ 1—\ + Kdl,[p,a) = — --+ Kdt,[p,a) 

< +Kdl{p,a) = (—±— + Kd,,ip,a))di{p,a), 



(o") " ' \4Ami„((T) 

where we used — cr||2 < — cr||i/ for the traceless matrix p — a. To upper bound the last term in 
parentheses by l/Amm(o'), we set e := 3 / {4K Xmini<^)) and the claim follows. □ 

Remark: A general linear upper bound of the form dsip, u) < Cdtr{p, a) as in Eqn. OJ cannot hold 
for cr ^ . For instance, in this case there exist density matrices p ^ Sd orthogonal to a (pa = 0). Then 
defining ps := 6p+ {1 — d)(T (for S e [0, 1]) one has dty{ps, a) — S and F{ps, <t) = — 5, so that 

dtr{ps,<7) < 2dl{ps,a) V(5e [0,1] . 

Thus, there do not exist constants C, e > such that dsip, a) < Cdti {p, cr) holds for all dtr(p, f ) < e. 



Appendix B: Proof of Theorem |3] 



This appendix proves Theorem[3]from Section HTCl 
Proof Let C be the matrix representation of the Liouvillian. The Jordan normal form gives 

S 0J,(A,) 



for some invertible matrix S, where A^ are the eigenvalues of C, and are Jordan blocks of the 

following form (note that eigenvalues Aj with Rc(Aj) = have one-dimensional Jordan blocks, so in 
particular all eigenvalues): 



A, A, 



A, A, 



Let dj > 1 be the dimension of Jordan block j. Then, in the Jordan basis {|fc)} defined by this. 



1=1 fe=i ^ ' 



(Bl) 



Let II • 1 1 be the operator norm, and denote by k{S) :=||S'||||S' ^|| the condition number of the similarity 
transformation into Jordan form. Then: 



S flft e"A^,)s-^ < K(5)e-*^ max 



dj I 



< K{S)e~'' max ( V V IMj.^' 



< Ce-*^max{(iA)"'-\l} 



\k){l\ 



where C is a t-independent constant and J := maxj dj is the dimension of the largest Jordan block. The 
last step is obtained by factoring (tX)''^^ out of the sum (for times t > 1/A) and bounding the remaining 
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term by a constant. Thus clearly, for any v <\ there exists a constant K > such that the last expression 
is upper bounded by A'e^*^ (for all t > 0). 

The lower bound is obtained similarly (letting Ji(Ai) be any Jordan block with Rc(Ai) — — A): 



T^W > Kisy 



e 



i:Rc(A,)#0 



k(S) ^ max 



M > i^iSy^ max 

M l<k<l<di 



,t./i(Ai 



kl 



-1 „-tA 



where maximum in the second-to-last expression runs over all matrix elements (fc, in dBll ). and in the last 
step we chose k ^ I ^ 1. 

Using these upper and lower bounds on ||Tt — 3^||, we invoke Lemma [12] (below) to complete the proof 
of Theorem[3] by lumping all of the time-independent constants together and denoting them by L > and 
R. □ 



Lemma 12 Let T : Add ^ be a quantum channel, and \ \ ■ \ \ the operator norm. Then: 



\T-%\\ < Vtr[T] < 



Vd 



\T-T^ 



(B2) 



Proof For the lower bound, we use in the first inequality below that every X € M.d can be written as 

X = Pi — P2 + iP-i — iPi with positive semidefinite Pi satisfying P1P2 = P3P4 = 0, and that then 
Il-P? ll2 < Il-P«ll2 = ll-'^lli- In the following chains we also use that | |X||2 < ||X||i < v^||X||2. 

||r-f^|| = sup ||(r-r^)(x)||2 < 4 sup ||(T-r^)(p)||2 

||X||2<1 P>0,||P||2<1 

<4 sup ||(r-r^)(P)||i = 4Vrf sup ||(T-T^)(P)||i = 8Vd77t,[T] . 

P>0,\\P\\l<Vd P>0,tr[P]<l 



1 

< - 
- 2 



sup ||(T- 
ll^lli<i 



1 

< - 
- 2 



sup 

Il-^ll2<l 



Vd\\{T^T^){X)\\, = ^\\f-f, 



□ 



